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ABSTRACT 

I derive up to second order in Eulcrian perturbation theory a new relation between the 
weakly nonlinear density and velocity fields. In the case of unsmoothed fields, density 
at a given point turns out to be a purely local function of the expansion (divergence) 
and shear of the velocity field. The relation depends on fi, strongly by the factor 
fiP) ~ g^^^ weakly by the factors K(Q) oc fi-^/es and C{n) oc f^-^/^i. The 
Gramann solution is found to be equivalent to the derived relation with the weak 
rj-dependence neglected. To make the relation applicable to the real world, I extend it 
for the case of smoothed fields. The resulting formula, when averaged over shear given 
divergence, reproduces up to second order the density-velocity divergence relation of 
Chodorowski & Lokas; however, it has smaller spread. It makes the formula a new 
attractive local estimator of large-scale density from velocity. 

Key words: cosmology: theory - galaxies: clustering - galaxies: formation - large- 
scale structure of the Universe 



1 INTRODUCTION 

The value of SI remains still one of the most intriguing un- 
knowns in cosmology today. The parameter Q,, defined as 
the ratio of the mean to the critical density is so crucial for 



cosmol q i g}' because its value deieimiiies the global geomeU y 
and ult imate fate of Univoroo. 



One way to measure S7 is to compare large-scale density 
fields of galaxies with the corresponding fields of galaxy ve- 
locities. Under widely accepted hypothesis of gravitational 
instability, observed large-scale peculiar flows of galaxies 
(deviations from Hubble flow) result from gravitational 
growth of initially small cosmic mass fiuctuations. The quan- 
titative relation between the mass density contrast field. 



5 = p/{(>) — 1, where {(>) is the tiieati density, atid the p6- 



culiar veluciLy field, i;, cau be deduced hum Lhe dynamical 
equations describing the pressureless self-gravitating cosmic 
fluid. For small density fluctuations linear theory can be ap- 
plied. Linear theory predicts that the density contrast is a 
linear and local function of the velocity divergence, 



(x). 



(1) 



In the above, the function f(Q) ~ fl'^'^ (see e.g. Pee- 



bles 19^ 0) and the superscript '(1)' denotes the linear theory 
limit. (Distances are measured here in km s~^, so the Hubble 
constant H = 1 in this system of units.) The above formula 
can be used to reconstruct from a large-scale velocity field 
the linear mass density field, up to an Sl-dependent multi- 
plicative factor f{Q). The comparison of the reconstructed 
mass field with the observed large-scale galaxy density field 



may therefore serve as a method for estimating fl and as a 
test for the gravitational instability hypothesis. 

Indeed, a strong correlation between the galaxy den- 
sity and v elocity divergence fields has been fou n d in ob- 
servation s (Dekel et al. 1993; Hudson et al. 1991:; Sigad et 
al. 1997). However, equation (|l|) assumes linear theory while 
the fields in question are weakly nonlinear. Smoothing of the 
fields, necessary among other things to reduce large individ- 
ual distance-estimation errors and the shot noise content, 
must be performed over a limited scale in order to optimize 
the information present in the finite-volume data. The po- 
tent algorithm for the mass density reconstruction from an 
observed radial velocity fleld currently e mploys a Gaussian 



smoothing length of 1000-1200 km s~^( pekel 1994 Pekel 
et al. 1997). At these scales, typical (rms) galaxy density 



fluctuations are of the order of several tens per cent, in con- 
tradiction with an underlying assumption of equation (|l]) 
that (5 <C 1. On the other hand, they are not in excess of 
unity, wherefore the name "weakly nonlinear" . The need for 
a weakly nonlinear generalization of linear formula (|l|) has 
been quickly recognized. The present potent algorithm uses 
the formula of Nusser et al. (1991), which is the Zel'dovich 
(1970) approximation expres sed in Eulerian coor d inates. 
However, N-body simulations ( Mancinelli et al. 1994; Sanon 



et al. 1997) have shown that though Zel'dovich approxima- 
tion does much better than linear theory equation (jl|), it 
still does not predict correctly the weakly nonlinear density- 
velocity relation (hereafter DVR). 

Weakly nonlinear regime is the regime of applicability of 
perturbation theory. To begin with, linear theory solutions 
for the density and velocity divergence fields that give rise to 
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linear equation (|l|) are nothing but perturbative series trun- 
cated at the lowest, i.e. first order terms. A natural way of 
extending linear DVR into weakly nonlinear regime is thus 
to take into account higher order perturbative contributions 
for density and velocity divergence. This has been done by 
Chodorowski & Lokas (1997a, hereafter CL97), who com- 
puted weakly nonlinear density-velocity divergence relation 
up to third order in perturbation theory. (Second order con- 
tributions were included already by Bernardeau 1992.) The 
resulting extension of the linear formula offers also a method 
for separating the effects of and possible bias between 
galaxy and mass distributions (CL97; Bernardeau, Chodor- 
owski & Lokas 1997, hereafter BCL). 

One might worry that the perturbative approximation 
to nonlinear DVR breaks down soon after the linear rela- 



tion d oes s o . Ifowever, N-body simul ations (BCL 



et al. ] 997 ; Chodorowski et al. 



1997) show the opposite: 
robust fit 



the perturbative formula is a very good robust fit to N- 
body results in the whole cosmologically interesting range 
of smoothing radii. 

Higher order perturbative solutions for density and ve- 
locity divergence are nonlocal. As a result, the relation be- 
tween weakly nonlinear density and velocity divergence at a 
given point is no longer deterministic. Still, since the spread 
comes exclusively from higher order contributions the two 
fields remain strongly correlated and the mean trend can 
serve as a useful local approximation to the true nonlocal 
DVR. This is exactly what has been calculated by CL97, 
who found that the formula for the conditional mean density 
given velocity divergence is given by the third-order polyno- 
mial in velocity divergence. The reverse case (mean velocity 
given density) has been calculated by Chodorowski & Lokas 
(1997b). W ork is in progress on the theoretical prediction for 
the spread ( Chodorowski, Lokas fc Polio 1997), as well as its 
measurement in N-body simulations (|BCL| ; |g hodorowski en 
al. 199^1). 

Summarizing the above in the statistical language, the 
perturbative polynomial in velocity divergence is an unbi- 
ased but non-zero variance local estimator of density. It is 
then natural to ask a question: among all unbiased local esti- 
mators of density from velocity, is it the minimum-variance 
one? It is quite unlikely. The fact that in weakly nonlin- 
ear regime density and velocity divergence at a given point 
are not related in a deterministic way does not exclude a 
possible existence of a purely local formula for density as a 



function of some derivatives of the velocity field. 



The 



irrotationality of the flow implies only that the tensor of ve- 
locity derivatives is symmetric, Vi^j = Vj^i. Hence, it has 6 
independent components, while the velocity divergence in- 
volves only 3 diagonal ones out of them. 

Indeed, already mentioned Zel'dovich approximation 
which does involve off-diagonal components is, like linear 
formula (|l|), deterministic. The formula based on Zel'dovich 
approximation is a biased estimator, but it is not based on 
rigorous Eulerian perturbation theory. The application of 
perturbation theory and the inclusion of the off-diagonal 
components in the velocity derivatives tensor may result in 
a local estimator which is both unbiased and has smaller 
variance than the relation of CL97. Why not simply zero- 
variance? To make the estimator applicable to the real 
world, smoothing of the fields has to be included. Smooth- 
ing is a nonlocal operation, hence it necessarily introduces 



some spread into any nonlinear DVR; I will discuss this in 
more detail later on. Still, since the density-velocity diver- 
gence relation is nonlocal already for unsmoothed fields, for 
smoothed fields it can be expected to have greater spread 
than the corresponding relation which is local when un- 
smoothed. 

This paper is devoted to constructing such a lower- 
variance estimator of density from velocity. In section 2 I 
derive up to second order in Eulerian perturbation theory 
a purely local relation between unsmoothed density and ve- 
locity fields. In section 3 I relate it to the Gramann (1993) 
solution. Subsequently, in section 4 I derive from it the un- 
smoothed density-velocity divergence relation. In section 5 
I generalize my relation for the case of smoothed fields. Sec- 
tion 6 is devoted to the comparison of the relation to N-body 
simulations. Finally, I summarize the results in section 7. 



2 DERIVATION 

Let us expand the density contrast in a perturbative series, 
+ . (2) 

In the above, 5^^^ denotes the p-th order perturbative contri- 



bution, w hich is of the order of (Fry 1984; Goroff et 



al. 1986). Introduce a variable proportional to the velocity 

(3) 



divergence 

d{x,t) = ~f{nyv^ ■v{x,t) 

and expand it as well. 



(2) 



q{3) 



(4) 



The linear theory solution, equation (|l|), gives 
5^^\x) = -ff'^^^x) . (5) 
From equations (^), (^ and (^) we have up to second order 
S{x) = ^{x) + (x) - {x) . (6) 

The second order perturbative contributions for S and 



■i? for ar bitrary S7 are (Bouchet et al. 1992 
al. 1995|) 



Bernardeau et 



and 



1 + K 



(5W)^ + V.5'i'.V.$W+. 



— $(^)<I.(^' (7) 
) ,1.1 ' 



(2) 



C{5'-''>f + V.<5^'^ ■ V.*^'' + (1 - C)^]-^$ 



(1) 



(8) 



Here, ^^^\x,t) is the linear gravitational potential satisfy- 
ing the Poisson equation 

A,$('' = 5^'' (9) 

and I use the Einstein summation convention. The weakly 
f2-dependent functions K and C are 



and 



(10) 



(11) 



The approximation ( |l^ ) is accurate to within 0.4 per cent 
in the range 0.05 < < 3 (|Bouchet et al. 1992|) and (0) to 
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within 2 per cent in the range 0.1 < < 10 (Bernardeau et 
al. 199^). 

Subtracting equation (j^ from (Q) we have 

S(2) _^(2) + 2^ (#^))^ '^ ^'^'-^'g -(12) 

Note that the terms V^S''^^ ■ Vaj^'^' have cancelled out. 
The well-known solution to equation (^) for the gravitational 
potential is 



dV 5'^''>{x') 
in \x — x'\ 



(13) 



Hence, the difference S^-^\x) — §^'^\x), in addition to the 
local term ~ contains a nonlocal term due to 

the gravitational potential, ~ 'Y^a^i^nS^-^^ {x')5''^\x") ~ 
Y,OL^L^'d^^\x')d'^^\x"). The right-hand side of equa- 
tion (|l^) is thus not a function, but a functional of i9. On 
the one hand, this means that up to second order in pertur- 
bation theory the velocity divergence field still contains full 
information about the density field. On another hand, how- 
ever, this information is distributed non-locally, that is the 
value of the density contrast at a given point x is only deter- 
mined by the values of the velocity divergence in all space. 
Shortly, the relation between the weakly nonlinear density 
and velocity divergence fields is deterministic, but nonlocal. 
Such a relation is of no use for cosmology, since we have 
at our disposal only surveys of limited volume, with uncom- 
plete sky coverage, uneven sampling, etc. The key point is to 
find a local relation between the density and velocity fields. 

One of the ways to do this is to approximate the 
nonlocal-in-i? term in equation (jl^ by its conditional mean, 
given the velocity divergence at a point x (Bernardeau 1992 



CL97). The resulting density-velocity divergence relation is 



then local by construction but clearly not deterministic; in 
other words it has some spread. The idea of the present pa- 
per is difi'erent: to recast the nonlocal-in-i? term to the form 
involving other local derivatives of the velocity field. 

Since the velocity field in weakly nonlinear regime re- 
mains irrotational it can be expressed as a gradient of some 
velocity potential <^r,, 



(14) 



Taking the divergence of the above equation and using equa- 
tion (111) shows that in linear regime the velocity potential is 
proportional to the gravitational potential. 



and equations (p^-(p^ yield 



(15) 



(16) 



Let us decompose the symmetric tensor of velocity 
derivatives into its trace (expansion), 6, and the traceless 
part (shear), Uij, 



where in general 
and 



V 



Vk,k ■ 



(17) 



(18) 



(19) 



The symbol Sij denotes the Kronecker's delta. Note by com 
paring with equation (^ that 



We have 

1 ^2 I 2 

where is the shear scalar 

2 



(20) 
(21) 
(22) 

From equations ( [l6| ) and (|^) taken at linear order we obtain 

^^-\r'(0'^'f + r'{<r'^'f- (23) 

Introducing the above equation to equation (^^ and us- 
ing (0) yields 

5(2)_^(2) ^ l±J£^l2£f-^iaW\^ 



o 

1 + K — 2C 2, (1)n2 
7, / ((^ ) ■ 



(24) 



Denote the rms fluctuation of the linear density fleld by 
e, = ((5*^')^). Perturbation theory predicts that 0^^'' ~ 
5(f) ~ so we have 9^ = (e'l^ + e'^^ + . . .f = {9^^^f + 
C(e^) and similarly for a^. Therefore, constructing weakly 
nonlinear density-velocity relation up to terms quadratic in 
e we can substitute {O^'^^f by 6^ and (a'^^)^ by a^. By doing 
so and introducing equation (^) into (^ we finally obtain 



S{x) = -f-^9ix)+ ^ + ^~^^ f-'e\x) 
- ii4T^rV(a.) + 0(e^). 



(25) 



This is the DVR computed up to second order in pertur- 
bation theory. The density is a local function of the expan- 
sion and the shear of the velocity fleld. A similar relation has 
been independently derived by Catelan et al. (1995), though 
without the decomposition of the term VijVij into the ex- 
pansion and the shear (cf. eq. [41] of Catelan et al. 1995). 
The relation depends on fl, strongly by the factor / and 
weakly by the factors K and C given by equations (p^-(pT[). 
Neglecting the weak f2-dependence we obtain 



5{x) = -f-'e{x) + ^r'e\x) - |r'<^' 



(26) 



Formula ( |25| ) for density involves only flrst derivatives of 
the velocity field. We owe this to the fortunate cancellation 
of the terms Vaj^'^"* • Vaj^'^-* in equation (p^. Had not these 
terms cancelled, we would have had on the right hand side 
of equation (M) an extra contribution ~ Va,"!?'-^) ■ Va;^'^'' ^ 
VjVi^ij -I- 0(e^7 involving the second derivatives. The practi- 
cal applicability of such a formula to extracting density from 
a very noisy velocity fleld would be doubtful. 



3 RELATION TO THE GRAMANN SOLUTION 

Gramann (1993) derived DVR up to second order in La- 
grangian perturbation theory. The form of this relation ex- 
pressed in Eulerian coordinates is 



5{x)^-f-^9{x) + -f-^m^{x)., 



(27) 
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where 

3 



(28) 



or, explicitly 

rriv — i;i,i'y2,2 + «i,i^^3,3 + «2,2«3,3 

— Wl,2?^2,l — Ul,3«3,l - U2,3«3,2 • (29) 

In the following I will relate equation ( p7| ) to my solu- 
tion for weakly nonlinear DVR calculated in the previous 
section. Let us decompose the tensor of velocity derivatives 
more generally into the expansion, shear, and vorticity (the 
asymmetric part), respectively. 



(30) 



1 

Here, tijk is the completely antisymmetric tensor; £123 = 
1. After applying the above decomposition, a few lines of 
algebra yield 



'^l.lW2,2 + '^l.ll'3,3 + W2,2l'3,3 



0"llC"22 



criiO"33 + cr220"33 
(31) 



2 2 1 2 

f^Vi + 0"23 ^ : 



Similarly, 

2 

f 1,2U2,1 + Wl,31'3,l + 'y2,3^'3,2 = Cri2 

where lj^ is the vorticity scalar, 
2 _ 

LU = Ijii ■ LJ — UdkL^k ■ 

By definition, 

a — (JijOij = 2 (cri2 + (713 + (723 ) + <^11+ 1722 

Using the above equation and the identity 

((711 + (722 + (733)^ = , 

we can recast equation (B2|) to the form 



2 

0-33 ■ 



(32) 

(33) 
(34) 
(35) 



(711(722 + (711(733 + (722(733 



I and (§ 



LJ . 

4 

we obtain 



Wl,2U2,l + Wl,3l'3,l + U2,3f3,2 

Using equations (^), ([ 

1 ^2 1 2 ,1 2 
3 2 4 

The final step is to use this result in equation ( p 
is 

5{x) = -r^e{x) + ±f-'9\x)-^f-'a\x) 



(36) 



+ 



1 r~2 I, \ 
Iff ^ W 



(37) 
, the result 



(38) 



Thus, the weakly nonlinear density in a given point is 
in general determined by the local values of the three scalars 
that can be constructed from the derivatives of the veloc- 
ity field: the expansion, shear and vorticity. Since before 
shell crossing a cosmic velocity field is irrotational we can 
in the above equation drop out the vorticity term. Then 
the equation coincides exactly with equation (^^ of the 
previous section. Equation ( Pq) is the DVR up to second 
order in perturbation theory with the weak ^-dependence 
neglected. Summarizing, the Gramann (1993) solution (|27|) 



is the second-order DVR with only the strong S7-dependence 
included. 

Equation ( p8[ ) bears some resemblance to the Ray- 
chaudhuri (1955) differential equation for the evolution of 
the velocity expansion. There, the source terms are simi- 
larly proportional to density, expansion, shear and vorticity 
(cf. eq. [22.14] of Peebles 1980). 



4 DENSITY — VELOCITY DIVERGENCE 
RELATION 

As already mentioned in section 2, density and velocity di- 
vergence at a given point are not related in a deterministic 
way. Let us rewrite equation (|25|) in the form 



S{x) 



'd{x) 



-i)\x) 



(39) 

and I have defined 
(40) 



where is related to 9 by equation 

Sij = —f o'ij . 

The spread in the S-iS relation clearly comes from the shear. 
The mean trend, defined as mean 5 given i9, is 

By definition of the conditional moment, 
_ / E2p(i9,E)dS 



(42) 



where p(i?, E) is the joint probability distribution function 
(PDF) for expansion and the shear scalar. It is sufficient to 
know the form of this PDF for linear 1? and E since already 



at the lowest order E^ 

-3\ 



e and nonlinear corrections are 



How to derive the joint distribution for i?'^' and E'^'? 
In the derivation of its general properties I will follow 
Juszkiewicz et al. (1995; Appendix A). I assume that the ini- 
tial conditions are Gaussian. Under this assumption, both 
t?*-^^ and five independent components of the shear tensor 
E^j^' are Gaussian distributed. Consequently, p(i9'-^', E^j^') 
is a multivariate Gaussian, entirely determined by its co- 
variance matrix. It is more convenient to compute the co- 
efficients of this matrix in the Fourier space. The Fourier 
transform of t?^^-' is obviously equal to the Fourier trans- 
form of the linear density field, ■i?^^'' = S'^^ . Thereafter I will 
drop out the superscripts '(1)'. The power spectrum P{k) is 
defined by the relation 



(43) 



where 5n denotes the Dirac's delta. The Fourier transform 
of a shear component is 



where ki = ki/k. We have 
d^fcd^fc' 



(27r)6 
_d^fe 
(^ 



/ - / - , 1 
iki kj — —i 



{SkSk')i 



(44) 



P{k) 



(45) 
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The second step uses equation (^3|) and the last one is obvi- 
ous by symmetry. It means that the linear shear components 
are uncorrelated with the linear velocity divergence. In gen- 
eral, the uncorrelation of random variables is only a neces- 
sary condition for their statistical independence; for Gaus- 
sian variables however it is also the sufficient one. There- 
fore, 1? and Sij are statistically independent; consequently 
the variables and E = (EijEi^)^''^ are also statistically 
independent. It implies 



p(^,E) =p(^)p(E). 

Using this result in equation (^2|) yields 

(E^>|^ = / E2p(E)dE, 



(46) 



(47) 



that is the mean value of E^ does not depend on ■!? and is 
equal to the ordinary mean. We have 

(E.,E,,) = ($,,$,,^)_i(52). (48) 

The integration by parts of the term yields 
{^,ii'^.n) = ('5^> = e^. hence 

(49) 

(cf. Silk 1974). Introducing this in equation ( ^ ) we finally 
obtain 



(<5)U = j9 + a2(^?' 



where 



a2 = 



1 + K -2C 



(50) 



(51) 



Equations (||)-(|ll) agree with the results of CL97 ob- 
tained by a completely different method, namely the Edge- 
worth expansion of the joint PDF for the variables 5 and 
■!?. Specifically, formula ( pl| ) for the coefficient 02 coincides 
with equation (81) of CL97 (for the case of no smoothing). 

I will now compute the spread around the mean trend, 
or the conditional variance. Hereafter in this section I will 
neglect the weak dependence of equations ( ^9| ) and (^ij) on 
fl. We then have 

((^ - wi-)')L = + ^ ■ (^2) 

It is known that E^ is x^-distributed with 5 degrees of 
freedom; the variance of the Gaussian variable underly- 
ing the distribution is 2/15e^ (Gorski 1988; cf. also Groth, 
Juszkiewicz & Ostriker 1989). Therefore 



= 2 • 5 ■ 



225 



8 4 
45" ■ 



(53) 



32 



hence finally 

{(S~{S)M')l = ^^e' + 0{e') . (54) 

The above formula exactly coincides with the result 
of Chodorowski, Lokas & Polio (1997; for the case of no 
smoothing), obtained by means of the Edgeworth expan- 
sion. 

And if initial conditions are non-Gaussian? In this case 
property (^) still holds as its derivation does not require 
any assumption about initial conditions. Therefore, and 
Eij remain uncorrelated. It however does not necessarily 



mean that they are statistically independent. Uncorrelated 
non-Gaussian variables may be, but equally well may not 
be, statistically independent; see Kendall & Stuart (1973) 
for the examples of both cases. Catelan et al. (1997) show 
that for a certain class of non-Gaussian models i9 and E are 
indeed independent. For this class of non-Gaussian mod- 
els the density-velocity divergence relation is therefore the 
same as for Gaussian initial conditions, equation (|5^). Only 
the spread around the mean trend is different since E^ is no 
longer x^-distributed. The form of the S--d relation for other 
non-Gaussian models remains to be investigated. 



5 EFFECTS OF SMOOTHING 

Ganon et al. (1997) test various approximations for weakly 
nonlinear DVR by the means of N-body simulations. Among 
the approximations considered is the Gramann (1993) so- 
lution. Since the Gramann solution is equivalent to equa- 
tion ( psj ) with the weak 57-dependence neglected (section 3), 
its properties essentially refiect the properties of the present 
solution. Ganon et al. (1997) plot the difference between 
the approximate and true density as a function of the true 
density, T> — 'D{5). The Gramann solution gives significant 
residuals which have a parabolic form. Is equation (051) thus 
incorrect? There is certainly no error in its derivation, but 
it cannot be straightforwardly applied to the smoothed den- 
sity and velocity fields that were estimated from N-body by 
Ganon et al. (1997). 

Inferring the fields from observations one has to intro- 
duce smoothing. The first reason for doing so is to reduce 
the effects of large individual distance-estimation errors and 
the shot noise. The second one is that only the field of ra- 
dial peculiar velocities is directly measured; to reconstruct 
from it the full three-dimensional velocity field we need to 
assume its potentiality. This assumption is valid only before 
trajectory crossing (where the Kelvin's circulation theorem 
holds), that is for the smoothed fields. However, even if we 
were provided by Nature with an accurate 3D velocity field 
(as accessible in N-body simulations) we would still prefer 
to smooth it, in order to reduce the nonlinearities. We would 
do so because linear or weakly nonlinear density and velocity 
fields are strongly correlated and the form of this correlation 
offers us a possibility to test the gravitational instability hy- 
pothesis and to measure the value of Q. 

Smoothing of the fields is realised by averaging them 
with a certain window of a certain scale R, Wr. For example, 
the smoothed density contrast field is 



5{x) 



A^x' 5{x')Wr{x 



(55) 



The velocity fields are commonly smoothed with a Gaussian 
filter. In this case, the second-order density-velocity diver- 
gence relation still has the form (^o|), but the coefficient 02 
is a function of the power spectrum index n, 



0,2 



1 + K ~2C 



•.Fi 



'n + 3 n + 3 5 1 
. 2 ' 2 ' 2' 4 



(56) 



Here, 2^1 is the hypergeometric function (CL97). The effec- 
tive power-law index n is the slope of the log P-log k relation 
at the smoothing scale R ( Bernardeau 1994| ). For n = —3 the 
coefficient 02 is {1 + K — 2C) /3, equal to that for the case of 
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no smoothing, equation (|5l|). The factor {l-\-K — 2C)/3 can 
be well approximated by its value for S7 = 1: 4/21 ~ 0.19. 
For higher spectral indices 02 grows monotonically up to the 
value ~ 0.30 for n = 1. At the scales of interest, i.e. of sev- 
eral megaparsecs, the effective index of the observed power 
spectrum is clearly different fr om the value w = — 3; e.g. for 
IRAS galaxies it is n = —1.4 (Fisher et al. 1993). Also, the 



power spectra which are commonly used in N-body simu- 
lations have the values of the effective index different from 
n — —3; e.g. for the standard CDM, for R > 5/i~^ Mpc, 
n > — 1 (e.g. CL97). Therefore, equation (|5c|), when applied 
to the smoothed fields, underestimates density because the 
value of the coefficient 02 is underestimated. Specifically, 



V = Aa2((5 



(57) 



where Aa2 = 12 (—3) — a2(n). The residual T> is thus 
a parabola in 5 with a negative coefficient. Finally, rela- 
tion (^Hj), equivalent to the Gramann solution, must yield 
the same residual, since relation ( |50| ) is its version av- 
eraged over the shear. This is indeed observed in the sim- 
ulations: the Gramann solution yields essentially the same 
parabolic residual as the S-i) approximation of Bernardeau 
(1992), derived also for the case of no smoothing. 

How to extend equation (p5|), or its another form (|39|), 
for the case of smoothed fields? Applying a smoothing filter 
to both sides of equation (pSj) we have 



5{x) 



= i9(a;) 
+ O(e'), 



1 + - 2C 



i92(a;) 



l + K ~2C 



(58) 
(59) 



However, from observations we can estimate only smoothed 
fields, and only after smoothing we can perform transforma- 
tions on them (like squaring). Our pur pos e is therefore to 
express the right-hand side of equation (Bq) as a function of 

-& and 



where e.g. the smoothed i)'^ is given by 
'¥{x)= d'x' '&^{x')Wr{x - x') . 



(60) 



Since smoothing and nonlinear transformations do not com- 



mute, in general i?^ jg j^ot equal to 1 
It is most clearly seen by writing 



= (i? - -df + / > 1?^ 



; similarly for the shear. 



(61) 



wherefrom i?^ jg equal to i9 only when the spectral index 
n — —3, i.e. when the fiuctuations have so large wavelengths 
that '& = '&. Moreover, the same values of i9 can lead to 
different values of (As a simplest academic example the 
reader can consider one-dimensional fields &i = 1 and &2 = 
2x, smoothed with a top-hat filter over the segment [0, 1].) 
It means that the relation between S{x) and the variables 
'&{x) and S(x) is non-deterministic. Again, the mean trend 
is given by the conditional mean. 



1+K 



2 



(62) 



The standard approach would be to derive the joint 
PDF for the four variables 1^, E^, ^, and E, p{W, E^, i9, S), 
and to integrate over it the second term on the right-hand 



side of the above equation. Fortunately, this horrible calcula- 
tion is unnecessary because the result can be simply guessed. 
Firstly, it must be a quadratic form in i9 and S, since it comes 
from second-order perturbative contributions. Second, when 
averaged over E, it must reduce to the second term of the 
smoothed version of equation (po|). 



-2\ 



where 02 is given by equation (|5 
I postulate that 



1 + K-2C ^^_3^\ 



and e = ( 5 



(63) 



(64) 



The first condition is obviously satisfied. Similarly to equa- 



tion (| 



we have 



d'^fc 
^2^ 



(4^ 



(65) 



Here, WR{k) is the Fourier transform of the window 
function. Thus, also the smoothed fields i9 and E are statis- 
tically independent. It implies that (E^)!;;^- = (E^) — (2/3)e^ 
(eq. tel), hence the second condition is satisfied as well. Fi- 
nally, an additional term cx i9 E on the right-hand-side of 
equation (^) would violate it, because (i?E)|;;j = i?(E) 7^ 0: 
the average of the shear scalar does not, unlike the average 
of its component Eij , vanish since E is positive-definite. 

Thus, the form on the right-hand-side of equation ( |64[ ) 
is the unique quadratic form in -d and E which, when aver- 
aged over E, reduces to the second term of equation (|63[). 
Therefore, postulated equation ( |64[ ) is indeed correct. Com- 
bining it with (k32) we obtain 



or, using equations ( po| ) and ( ]4C| ) 

-2 /772 



3_2 

r 



(66) 



(67) 



Here, 9 and a is the expansion and shear of the smoothed 
velocity field v, 



and 

_ _ /_ _ -.1/2 



with 



(68) 



(69) 



(70) 



Equation (p6[), with the coefficient 02 given by equa- 
tion (pq), constitutes an extension of equation ( p9[ ) for the 
case of smoothed fields. As already discussed, smoothing 
of the fields induces spread in the relation between the 
smoothed density and the smoothed expansion and shear 
in a given point. In the present paper I will not attempt to 
compute the spread. However, it is certainly smaller than 
the spread in the smoothed density-velocity divergence re- 
lation (|63|), since this relation is obtained from (^^ by aver- 
aging over the shear. This averaging is an extra source of the 
spread in the S-i!) relation: for unsmoothed fields the spread 
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is given by equation (p^, while the relation ( p9| ) between 
density and the two velocity scalars is entirely determinis- 
tic. 

It should be stressed that equation (^6|) assumes Gaus- 
sian initial conditions, since the coefficient a2, equation (pij), 
has been computed by CL97 under this assumption. For non- 
Gaussian initial conditions, a detailed form of the relation 
between the smoothed density and the smoothed expansion 
and shear remains to be studied. 



6 COMPARISON TO N-BODY SIMULATIONS 

Ganon et al. (1997) have performed N-body simulations for a 
CDM family of models. For the fields smoothed with a Gaus- 
sian window of radius R = 12h~^ Mpc, Ganon et al. have 
estimated the value of the coefficient 02 to be 



(NB) 



0.28 



(71) 



(cf. also Willick et al. 1997). As explained in the previous 
section, the Gramann approximation yields 



AG) 



21 



tT2r 



What is the value of 02 predicted by the approxima- 
tion derived here? The hypergeometric function in expres- 
sion (^^ can be expanded in powers of n + 3, simply by 
rearranging the appropriate terms of the hypergeometric se- 
ries (cf. Lokas et al. 1995). In the expansion, the term linear 
in n 4- 3 vanishes. The resulting formula for the coefficient 
a2, truncated at the cubic term, is 



a2 



= — [1 + C2(n + 3)^ + C3(n + sf] , 



where 



C2 



0.02139 



and 



C3 = 0.00370 . 



(73) 



(74) 



The polynomial ([r3|) approximates expression (jS^) in the 
range —3 < n < 1 with accuracy of 0.2% or better. Ifflj^l, 
the factor 4/21 should be replaced by the factor (1 + A' — 
2C)/3, weakly varying with Q. 

For a smoothing radius of 12h~^ Mpc, the effective in- 
dex of the standard CDM spectrum is (for details see CL97) 

n = -0.404. (75) 

Using formula (^), the predicted value of the coefficient 02 
is then 



a2 



0.23. 



(76) 



This value lies half-way between Oj*^' of the Gramann ap- 
proximation, equation (^2|), and a^^^^ estimated from N- 
body by Ganon et al. (1997), equation (|7l|). It means, on the 
one hand, that the density-velocity relation derived in this 
paper, equation ( p7| ) with 02 given by equation ( |56| ) or ([73|), 
is not in full agreement with the results of N-body by Ganon 
et al. On the plot of the difference between the approximate 
and the true density as a function of the true density, in- 
ferred from N-body, it will still leave some parabolic resid- 
uals. On the other hand, however, these residuals will be 
smaller, roughly two times, than for the Gramann approx- 
imation. As stated earlier, the Gramann approximation is 
valid only for unsmoothed fields. Not surprisingly then, the 



proper second-order formula for smoothed fields is a better 
estimator of smoothed density from smoothed velocity. 

The slight discrepancy between the value ( |7^ ) and that 
estimated from N-body can be attributed to perturbative 
contributions of higher orders. Here, or in the paper of CL97, 
the value of a2 has been derived at the lowest relevant order 
in perturbation theory, the second. Higher-order corrections 
yield a contribution to 02 which is of the order of = (5 ). 
In other words. 



a2(e) = a2(e ^ 0) + pe' + ©(e"* 



(77) 



where the dimensionless coefficient p remains to be calcu- 
lated. Since the (linear) variance of the field in question 
(Gaussian smoothed with radius of 12h~^ Mpc) is = 
0.076, the value p ~ 0.7 is sufficient to account for the dis- 
crepancy between a2 and Oj'^^^ 

On the other hand, it is important to confirm the results 
of Ganon et al. (1997) by independent simulations. Prelim- 
inary results of such show that the coefficient 02 has indeed 
slightly higher value than predicted by s econd order pertur 



bation theo ry, approaching it for small e ( BCL 



et al. 1997) 



Chodorowski 



7 SUMMARY 

In the present paper I have studied the relation between the 
weakly nonlinear cosmic density and velocity fields. I have 
derived up to second order in perturbation theory an expres- 
sion for density as a purely local function of the expansion 
(divergence) and shear of the velocity field (eq. |^ or ps]]). 
The relation depends on Q both strongly and weakly. I have 
shown that the Gramann (1993) solution (eq. ^^) is equiva- 
lent to equation ( |25| ) with the weak ^-dependence neglected. 
Also, a similar relation has been independently found by 
Catelan et al. (1995). 

The locality of the relation derived here is in contrast 
with the non-locality of the density-velocity divergence re- 
lation calculated by CL97. I have shown that averaging of 
equation (^^ over shear given divergence yields up to sec- 
ond order the formula of CL97 (eq. ISCl]); thus, the source of 
the spread in the latter relation is the distribution of shear. 
Since the form of this distribution is known, it is straight- 
forward to compute the higher-order conditional moments 
of the density-velocity divergence relation. In particular, I 
have computed the spread of this relation (the conditional 
variance; eq. |Q) and have found it to coincide with the 
result of rather cumbersome calculations by Chodorowski, 
Lokas & Polio (1997). 

Smoothing is a necessary ingredient of large-scale 
density-velocity comparisons. I have then generalized equa- 
tion ([^) for the case of smoothed fields (eq. j67| with 02 
given by eq. |H^). Smoothing not only modifies the shape 
of the relation between density and the two velocity scalars 
but also induces the spread in it. I have not computed the 
spread explicitly. Still, I have shown that it is smaller than 
the corresponding spread in the density-velocity divergence 
relation, since that spread has an extra source: averaging 
over shear given divergence. 

Equation (|2^) for unsmoothed fields does not depend 
on the type of initial conditions. Its smoothed counter- 
part (^^, however, assumes Gaussian initial conditions; a 
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detailed form of it for non-Gaussian models remains to be 
investigated. 

Checks against N-body simulations show that rela- 
tion (H) is still a slightly biased estimator of large-scale 
density from velocity (section 6). However, it was not the 
main concern of this paper. The ultimate goal of attempts 
to establish semi-linear relations between the density and 
velocity fields is to construct a local estimator of density 
from velocity which is not only unbiased but has minimum 
variance as well. The present paper is the first attempt to 
address the second point. As stated earlier, the inclusion of 
the shear term in the relation between the large-scale den- 
sity and velocity reduces its spread. In this paper I have 
demonstrated how to include the shear in the second-order 
perturbative relation. In future, I plan to do so for higher- 
order relations as well. 

The relation density versus velocity expansion and 
shear, derived here, is a new, lower-variance estimator of 
density from the velocity field that can be applied in the 
cosmic density-velocity comparisons. 
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